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Abstract. Discrete geometric estimators aim at estimating geometric 
characteristics of a shape with only its digitization as input data. Such 
an estimator is multigrid convergent when its estimates tend toward the 
geometric characteristics of the shape as the digitization step h tends 
toward 0. This paper studies the multigrid convergence of tangent esti- 
mators based on maximal digital straight segment recognition. We show 
that such estimators are multigrid convergent for some family of con- 
vex shapes and that their speed of convergence is on average ©(/is). 
Experiments confirm this result and suggest that the bound is tight. 

1 Introduction 

The problem of estimating geometric quantities of digitized shapes like area, 
perimeter, tangents or curvatures is an important and active research field. Dis- 
crete geometric estimators have indeed many applications in vision, shape anal- 
ysis and pattern recognition. It is however difficult to compare objectively their 
respective accuracy, since for a given shape there exists infinitely many shapes 
with the same digitization. In this paper, we are mainly interested by the multi- 
grid convergence property of some estimators, which is one of the few existing 
objective criteria. This property ensures that a better resolution brings a better 
approximation. For concrete applications, the speed of convergence is a very im- 
portant criterion too, since it has a huge impact on their accuracy at standard 
resolutions. 

Formally, taking the definitions of [9], let Dig^ be some digitization process 
of step h. Let F be a family of shapes in and let G be a geometric quantity 
defined for all X €: V. A discrete estimator of G is a map that associates to a 
digitization Digf^{X) an estimation of G{X). A discrete estimator Q is multigrid 
convergent toward G for F and Dig iff, for any A" 6 F, there exists some hx > 
for which 

yO<h< hx, |e(Dig,(A)) - G(A)| < T{h) , 

where t : R+ — > R+* has limit value at ft, = 0. This function defines the speed 
of convergence of G toward G. 

For instance, for the family of plane convex bodies with C'^-boundary and 
positive curvatures, denoted later on by F^, the best known speed of convergence 
for an area and or a moment estimator is 0(/itt) |8ll0j . for a perimeter estimator 
it is 0{h) [l6]. There are fewer results concerning local geometric quantities like 



tangent direction or curvature. The first works on this topic were presented in [3], 
where some evidence of the convergence of tangent estimators based on digital 
straight segment (DSS) recognition were given. We report the recent work of 
[13] which estabhshes the muhigrid convergence of tangent direction estimators 
based on maximal DSS recognition, with an average speed of convergence of 
0{h^). The recent result of ^ has confirmed that there is yet no curvature 
estimator proven to be multigrid convergent. 

In this paper, we prove a new upper bound for the average speed of conver- 
gence of discrete tangent estimators which are based on maximal DSS recognition 
around the point of interest |6ll2j . This new bound of ©(/is) was suggested by 
the experimental study of [T^] and enhances the previous bound of ©(ft 3) [TJ]. 
The proof of this enhanced result, obtained for shapes in F^, follows these steps: 

Section 2 Digitizations of convex shapes are digital convex polygons (CDP). 
We achieve thus a better localization of the shape boundary with respect to 
the digitized boundary (Proposition [1]). 

Section 3 We first claim that the DSS characteristics have an average asymp- 
totic behaviour depending on their length (Claim [J). The asymptotic edge 
length of CDP is recalled (Theorem [T]), which induces a superlinear localiza- 
tion of the shape boundary in 0{h^) (Proposition [2]). With these relations, 
the slope of digital edges of CDP is shown to be multigrid convergent to 
the tangent direction of a nearby boundary point with an average speed of 
0{hi) (Proposition El). 

Section 4 The behaviour of maximal DSS is identical to the behavior of digital 
edges, since any maximal DSS contains at least one digital edge (Proposi- 
tion 2]). The average speed of convergence of tangent estimators based on 
maximal DSS is thus achieved (Theorem [3]) . 

We will then conclude and open some perspectives in Section 5. 

2 Preliminary definitions and first properties 
2.1 Digitization and convex digital polygon 

Let S be some subset of R^. Its Gauss digitization of grid step h is defined as 
D^(5) = S* n /iZ X KL. Thus, the considered digitized objects are subsets of the 
rescaled digital plane /iZ x hi. 

A convex digital polygon (CDP) P is a subset of the digital plane hT, x KL 
with a single 4-connected component equal to the Gauss digitization of its convex 
hull, i.e. r — Dj^(conv(/^)). Its vertices (l^)i=i..e form the minimal subset for 
which r — D^(conv(Vi, . . . , Ve)) (vertices are ordered clockwise). The points of 
r which are 8-adjacent to some point not in F form the border of F. It is a 
4-connected digital path that visits every Vi in order. When moving clockwise, 
any subpath /Vi,Vi+i is called a digital edge of F, while the Euclidean straight 
segment ViVi+i is an edge of F. 

For small enough grid steps. Gauss digitizations of any shape in are convex 
digital polygons. This result will be admitted throughout the paper. 



2.2 Standard line, digital straight segment, maximal segments 

Definition 1. JlSf The set of points {x,y) of the digital plane J? verifying \i < 
ax — by < fi + \a\ + \b\, with a, b and jjL integer numbers, is called the standard 
line with slope a/b and shift fi. Standard lines are the ^-connected discrete lines. 

The quantity ax — by is called the remainder of the line. Points whose re- 
mainder is fi (resp. |a| + |6| — 1) are called upper (resp. lower) leaning points. 
Any finite connected portion of a standard line is called a digital straight segment 
(DSS). Its characteristics are the slope a/b and the shift fi of the standard line 
containing it with smallest \a\ + \b\. 

Most of the results demonstrated here are directly transferable for 8-connected 
curves since there is a natural bijective transformation between standard and 
naive digital lines. In the paper, all the reasoning is made in the first octant, but 
it extends naturally to the whole digital plane. 



2.3 Recursive decomposition of DSS 

We here recall a few properties about patterns composing DSS and their close 
relations with continued fractions. They constitute a powerful tool to describe 
discrete lines with rational slopes [2|7| . Without loss of generality all definitions 
and propositions stated below hold for standard lines and DSS with slopes in 
the first octant (e.g. |- with < a < 6). In the first octant, only two Freeman 
moves are possible, is a step to the right and 1 is a step up, 4-connected digital 
paths can be expressed as words of {0, 1}*. 

Definition 2. Given a standard line (a,b,fj,), we call pattern of characteris- 
tics (a, b) the word that is the succession of Freeman moves between any two 
consecutive upper leaning points. The Freeman moves defined between any two 
consecutive lower leaning points is the previous word read from back to front and 
is called the reversed pattern. 

Since a DSS has at least either two upper or two lower leaning points, a DSS 
(a, 6, ^) contains at least one pattern or one reversed pattern of characteristics 
(a, 6). There exists recursive transformations for computing the pattern of a 
standard line from the simple continued fraction of its slope (see [2] , |18j Chap. 4 
and [S] Chap. 9), here Berstel approach better suits our purpose. 

A rational slope z in ]0, 1] can be written uniquely as the continued fraction 
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and is conveniently denoted [0, ui, . . . ,u„]. The Ui are called the partial coeffi- 
cients and the continued fraction formed with the k first partial coefficient is 
said to be a k-th convergent of z and is a rational numbers denoted by Zk- The 
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Fig. 1. A digital straight segment of characteristics (10, 23, 0) with an odd depth 
slope, taken between the origin and its second lower leaning point 



depth of a k-th convergent equals fc. We conveniently denote pk the numerator 
(resp. qk the denominator) of a k-th convergent. 

We recall a few more relations regarding the way convergents are related and 
which will be used later on in this paper: 

Vfc > 1 Vk<lk-x - Pk-iQk = (1) 
Pq = p-i = l Vfc > 1 Pk = UkPk-1 + Pk-2 (2) 

go = 1 <?-i Vfc > 1 qk= UkQk-i + qk-2 (3) 

Given a rational slope between and 1 its continued fraction is finite and for 
each i, m is a strictly positive integer. In order to have a unique writing we 
consider that the last partial coefficient is greater or equal to two; except for 
slope 1 = [0, 1]. 

Let us now explain how to compute the pattern associated with a rational 
slope z in the first octant. Let us define E a mapping from the set of positive ra- 
tional number smaller than one onto the Freeman-move's words. More precisely: 
E(zq) — 0, E{zi) — 0"il and the other values are expressed recursively: 

E{Z2^+1) - E{z2^r^' + 'E{z2^^l) (4) 
E{Z2^) = E{z2^-2)E{Z2^-Ir^' (5) 

It has been shown that this mapping constructs the pattern (a, b) for any 
rational slope z = |. Fig. [T] illustrates the construction of an odd pattern using 
the mapping E. The Minkowski length of E{zk) equals pk -\- qk and can be 
computed recursively using Eq. ^ and ([3]). Moreover we recall that any digital 
edge is a pattern or a succession of the same pattern, its digital parameters are 
its slope, denoted by |, and the number 5 of repetitions of the pattern E{^). 

2.4 Localization accuracy of digitized convex shapes 

One can expect that the boundary of a convex shape is approximately at distance 
h from the border of its digitization of grid step h. In fact, for a convex shape 



Fig. 2. Two local constraints for the real underlying convex shape 



5, its boundary dS is much closer to the convex hull of 0^(5) than h. A better 
localization of the shape is thus possible, as stated below and illustrated in Fig. [21 

Proposition 1. Let S be a convex shape such that D^(S') is a CDP F for some 
h. Consider an edge V^Vi+i of F with slope in the first octant. Then any point of 
the boundary dS above the straight segment ViVi+i has a vertical distance to it 
no greater than ^ ^ ^ , where the slope of this edge is the irreducible fraction ^ . 

The proof of this proposition can be found in 

3 Asymptotic behaviour of edges of digitized shapes 

We study here asymptotic properties of edges of digitized convex shapes. Their 
average length is first exhibited and an experimental study of the digital param- 
eters of edges is presented. The direction of digital edges is shown to converge 
toward the tangent direction. 

3.1 Asymptotic number and digital parameters of edges 

Let S be some shape of Fj?. We have the following theorem: 

Theorem 1. (Adapted from [JJ, Theorem 2) For a small enough h, the Gauss 
digitization 0^(5*) of S is a CDP and its number of edges ne(D^(S')) satisfies: 

ci(5)-V<n,(DG(5))<C2(5)-^ 

ha h'i 

where the constants ci{S) and C2{S) depend on extremal bounds of the curvatures 
along S . Hence for a disc Ci and C2 are absolute constants. 

As an immediate corollary, the average Minkowski length of edges grows 
as 0{ha). The question is: what is the average behavior of the digital parameters 
p, q and 6 of edges, knowing that the average digital length I = S{p + q) satisfies 
0{hi/h) = 0{l/m) ? Since the slope of edges should tend toward the slope of 
points on the boundary of S and since almost all these points have irrational 
slope, p and q should tend toward infinity almost everywhere with a bounded 
5. This observation is confirmed by experiments, as illustrated on Fig. [31 which 
plots the means and standard deviations of j and 5 for edges on finer and finer 
digitizations of a disk. It is clear that q (but also p) satisfies the same asymptotic 
law as / while 5 remains bounded on average. We make hence the following claim: 
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Fig. 3. Plots in log-space of the means of digital parameters j (left) and S (right) 
for the edges of Gauss digitizations of a disk, as functions of the inverse i of the 
grid step. Standard deviations are symbolized with crrorbars 



Claim 2. On average, the digital parameters p, q of the edges o/D^(S') with 
5* G grow as 0{^), whereas 6 is bounded, when h tends toward 0. 

This claim induces a first result about the convergence speed of the localiza- 
tion constraints of Proposition [T] {qn-i has the same asymptotic law as q„). 

Proposition 2. Let S E and D^{S) its Gauss digitization. For a small 
enough h, in each octant, the vertical distance between any edge ofD'f^{S) and 
dS is bounded. On average, this bound is 0{h^) (and never worse than 0{h)). 

3.2 Convergence of tangent estimation based on edge direction 

Let M be a point on the boundary of a shape 5 G and let 9m{S) be the 
tangent direction of dS at AI. We propose to estimate d]\i{S) as the arctangent 
of the slope of the digital edge of the CDP 0^(5") lying below M (for the first 
octant; definitions for other octants are straightforward). We denote this discrete 
estimator by T^f^'^. Assuming Claim[2l we have the multigrid convergence of this 
tangent estimator: 

Proposition 3. Let M be a point in the plane and let F^(Af) be the subset 
of shapes S of F^ with M G dS . The discrete estimator T^^ ^'^ at point M is 
multigrid convergent toward the tangent direction Om for¥^{M) and Gauss dig- 
itization. Furthermore, its speed of convergence is on average 

Proof. We have to prove that for a shape S G Fj?(Af), there exists some grid 
step hs for which VO < /i< hs, \Tlf^''{V)f,S) - 6'a/(S')| < t(/i). Without loss of 
generality, we assume that the tangent direction 9m{S) is in the first octant and 
we locally parameterize dS as (a;, f{x)), setting M as the origin of the coordinate 
axes {{xm^Vm) = (0,0)). Let Hq be the digitization step below which 0^(5") is 
a CDP. 



Since the slope of dS around M is in the first octant, there exists some 
hs < ho such that for any < h < hs, the vertical straight segment going 
down from M intersect some edge ViVi+i of D^{S). Let 6, p and q be the digital 
parameters of this edge, let z be its slope - and let n be the depth of z (i.e. 
z — The horizontal distance between M and either Vi or V^+i is necessarily 
greater than half the horizontal length of ViVi+i. We assume Vi+i maximizes 
this distance, without loss of generality. We get: 

% < x.,^^, < Mqn (6) 

From Proposition [1] we have for any x £ [xy-,Xy-^-^^]: 

zx^^<fix)<zx+^ . (7) 

Inserting Taylor expansion of f{x) about a; = in Eq. ([7]) induces: 

\z-fiO)\<^ + Oix). (8) 

Setting X = x^i^i in Eq. ^ and using both sides of Eq. ^ gives the inequahty 

k-/'(0)|<7j3^ + O(Mg„) . (9) 

We notice that z = tan(7^'^^''(D^(S'))) and that /'(O) = tan(6'M(5')). Moreover, 
for any angle s,t £ [0, j], we have |s — i| < | tan(s) — tan(t)|. With these two 
remarks, Eq. ([9]) implies: 

I^M '°(D^(^)) - ^A/(^)l < id::^ + oihSq^) . (10) 

Since Sqn is no smaller than half the edge length, it follows from Theorem [1] that 
the dominant term j^— | — - is at least some on average. In jl3j it is also 

shown that there is no edge of bounded length as h tends toward 0. Then the 
right part of Eq. (jlOp tends toward 0. The multigrid convergence is thus shown. 
At last, assuming Claim [2l Eq. (fTO| then induces 

\T^^"'iBf,iS))-eMiS)\<Oihi) , (11) 

which indicates that the average speed of convergence of r^^s° is 0(/ii). □ 

4 Tangent estimators based on maximal DSS recognition 

This section discusses the convergence speed of discrete tangent estimators based 
on maximal digital straight segment (maximal segment) recognition. Along a 
digital path, maximal segments are the inextensible digital straight segments, 
otherwise said adding the next point to the front or to the back constructs a set 
of digital points that no standard line contains. The set of all maximal segments 
of a digital path can be extracted efficiently in time linear with its number of 
points [6] . Maximal segments have deep links with edges of convex hulls [5|14j . 




Fig. 4. For both plots, the digitized shape is a disk of radius 1 and the abscissa 
is the inverse of the digitization step. Left: plot in log-space of the mean absolute 
error between the A-MST tangent direction and the theoretical one (mean taken 
over all digital points). The convergence speed on this shape is likely to be in 
0{h^). Right: plot in log-space of the digital length of maximal segments, which 
is on average (9(l//i3) 



Estimating the tangent direction at some point is then achieved by consid- 
ering specific DSS [IT] or maximal segments [6] containing this point. A recent 
experimental evaluation [12] has shown that tangent estimators based on maxi- 
mal segments are accurate and preserve convexity properties of the real shape. 

Let M be a point on the boundary of a shape S* S and let 9m{S) be the 
tangent direction of dS at M . We propose to estimate 6m{S) as the arctangent of 
the slope of any maximal segment of the CDP D^(S') lying below M (for the first 
octant; definitions for other octants are straightforward). We denote this discrete 
estimator by T^'^. We shall prove in Theorem [3] that this estimator is multigrid 
convergent with average convergence speed oiO{hi). As a corollary, the Feschet- 
Tougne tangent estimator [6], which uses the most centered maximal segment 
around Af , and the A-MST estimator [12], which makes a convex combination of 
the directions of all maximal segments around Af, are also multigrid convergent 
with same average speed. 

Before proving this theorem, note first that the experimental evaluation of 
the A-MST estimator, whose absolute error is displayed in Fig.jH left, suggests 
that its average convergence speed is indeed likely to be in O^k's). Secondly, the 
average Minkowski length of maximal segments behave as the average length of 
the edges, that is in ©(/la) as shown on Fig. [H right (the Minkowski length is 
h times the digital length). Thirdly, every maximal segment contains a digital 
edge as stated in the following proposition whose proof can be found in [llj : 

Proposition 4. On the border of a CDP F , any maximal segment, whose slope 
is some contains at least one digital edge that has either the same slope or 
its n — 1- convergent . 

We may now state our main result, which relies also on Claim |2j 



Theorem 3. Let M be a point in the plane and let Fj?(A/) be the subset of 
shapes S o/Fj? with M G dS . The discrete estimator at point M is multigrid 
convergent toward the tangent direction 6m for Fi?(A/) and Gauss digitization. 
Furthermore, its speed of convergence is on average 0{h's). 

Proof. We take the same notations as in the proof of Proposition[3l We consider 
a maximal segment MS below M with slope z„. According to Proposition [H it 
contains a digital edge VV of slope z„ or Zn-i, which may not be below M. 
Worst case is for a digital edge with slope z„_i and this is the one considered 
below. We have: 

< |2n-l-/'(0)| + |z„-Z„_iI 

< |z„_i_/'(x)|+0(a;) + ^-^ , 

using Taylor relation and Eq. ^ to get the last inequality. 

If we choose some abscissa x such that (x, f{x)) is above the edge VV, then 
the proof of Proposition[3]indicates that \zn-i—f'{x)\ < 0{h3) on average. Since 
{x, f{x)) is above MS too, x cannot be greater than the length of MS, which is 
also some on average. At last, Claim [2] similarly provides — ^ — = 0{hl). 

Summing all these bounds concludes the proof. □ 

As an immediate corollary, the A-MST estimator [12] and the Feschet-Tougne 
tangent estimator |6' have the same asymptotic behaviour. 

5 Conclusion 

We have studied properties of digital convex polygons and exhibited several new 
results (Proposition[T]and Proposition|3]) . We have also examined the asymptotic 
properties of digital edges on digitized convex shapes (Claim[2]), which has led to 
a position estimator of average convergence speed ©(/i 3 ) (Proposition[2]) and to a 
tangent direction estimator of average convergence speed 0(^13 ) (Proposition[3]). 
At last we have shown the new bound of 0{h^) for the average convergence speed 
of tangent estimators based on maximal segments (Theorem which matches 
experimental evaluation. 

These results indicate that curvature estimators relying on surrounding DSS 
have most probably an error of 0{h^): position = f{x) in 0{h'^), tangent di- 
rection = f'(x) in 0{h3), curvature = f"{x) probably in 0(1) (uncertainty 
on tangent further divided by ©(/la)). The problem of exhibiting a multigrid 
convergent curvature estimator is thus still open. 

Another straightforward extension of this work would be to investigate the 
properties of discrete surfaces and estimators based on digital plane recognition. 
However since the problem of finding an enclosing polyhedron with a minimal 
number of 2-facets has been proven to be NP-hard (see [3]), the problem would 
get much harder than the two dimensional case. 
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